Conventional Forces can Explain the Anomalous Acceleration of Pioneer 10 
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Anderson, et al. find the measured trajectories of Pioneer 10 and 11 spacecraft deviate from the 
trajectories computed from known forces acting on them. This unmodelled acceleration (and the 
less well known, but similar, unmodelled torque) can be accounted for by non-isotropic radiation 
of spacecraft heat. Various forms of non-isotropic radiation were proposed by Katz, Murphy, and 
Scheffer, but Anderson, et al. felt that none of these could explain the observed effect. This paper 
calculates the known effects in more detail and considers new sources of radiation, all based on 
spacecraft construction. These effects are then modelled over the duration of the experiment. The 
model reproduces the acceleration from its appearance at a heliocentric distance of 5 AU to the last 
measurement at 71 AU to within 10 percent. However, it predicts a larger decrease in acceleration 
between intervals I and III of the Pioneer 10 observations than is observed. This is a 2 sigma 
discrepancy from the average of the three analyses (SIGMA, CHASMP, and Markwardt). A more 
complex (but more speculative) model provides a somewhat better fit. Radiation forces can also 
plausibly explain the previously unmodelled torques, including the spindown of Pioneer 10 that is 
directly proportional to spacecraft bus heat, and the slow but constant spin- up of Pioneer 11. In 
any case, by accounting for the bulk of the acceleration, the proposed mechanism makes it much 
more likely that the entire effect can be explained without the need for new physics. 

PACS numbers: 04.80.-y,95.10.Eg,95.55.Pe 



I. INTRODUCTION 

In [1], Anderson et al. compare the measured trajec- 
tory of spacecraft against the theoretical trajectory com- 
puted from known forces acting on the spacecraft. They 
find a small but significant discrepancy, referred to as 
the unmodelled or anomalous acceleration. It has an 
approximate magnitude of 8 x 10~* cm s^^ directed 
approximately towards the Sun. Needless to say, any 
acceleration of any object that cannot be explained by 
conventional physics is of considerable interest. These 
spacecraft have been tracked very accurately over a pe- 
riod of many years, so the data are quite reliable, and 
the analysis, though complex, has been reproduced by 
Markwardt [2] . Explanations for the acceleration fall into 
two general categories - either new physics is needed or 
some conventional force has been overlooked. 

One of the most likely candidates for the anomalous 
acceleration is non-isotropic radiation of spacecraft heat. 
This is an appealing explanation since the spacecraft dis- 
sipates about 2000 watts total; if only 58 watts of this to- 
tal power was directed away from the sun it could account 
for the acceleration. The bulk of the spacecraft heat is 
radiated from the two Radioisotope Thermoelectric Gen- 
erators (RTGs), which convert the heat of decaying plu- 
tonium to electrical power to run the spacecraft. The 
remainder of the heat is radiated from various spacecraft 
components as a result of electrical power dissipation, 
and by a few small Radioisotope Heater Units (RHUs) 



'Electronic address: louOcadence . com 



which serve to keep crucial components warm. At least 
three mechanisms been proposed that could convert heat 
radiation to net thrust - non-isotropic radiation from the 
RTGs themselves, heat from the RTGs reflected off the 
antenna, and non-isotropic radiation from the spacecraft 
bus. Anderson et al. reply with arguments against each 
of the proposed mechanisms. 

Although less well known, Anderson et al. also report 
that both Pioneers experience anomalous angular accel- 
erations. Pioneer 10 is spinning down at a rate corre- 
sponding to a torque of approximately 4.3 x 10~^ n — 
meters (in 1986). This torque is slowly decreasing - for 
most of the data span, intervals I and III of Anderson, 
the torque is directly proportional to the power dissipated 
by the spacecraft bus (in interval II, it appears that gas 
leaks dominate the spin behavior). This proportionality, 
and the size of the effect, lead naturally to an explanation 
of non-symmetric radiation of bus heat, supplemented by 
somewhat larger gas leaks in interval II. Pioneer 11, when 
not maneuvering, was slowly and constantly spinning up. 
The authors speculate that the source could be gas leaks. 

This paper argues once again that non-isotropic radia- 
tion is the most likely cause for both the unmodelled ac- 
celeration and at least some of the unmodelled torques. 
Each of the radiation asymmetries is re-examined, and 
a few previously unmodelled forces are included. Their 
sum is more than enough to account for the acceleration, 
and provides a plausible explanation for the unmodelled 
torques. Furthermore, we compare the acceleration in- 
duced by the proposed mechanisms with the measured 
data. We get reasonable, but not perfect, agreement over 
the whole data span. The main discrepancy is that the 
radiation thrust is predicted to decrease more quickly 



than the observed acceleration. The discrepancy is small 
(less than Icr) from the analysis of Markwardt [2] , but 
roughly a 2a discrepancy from the average results of the 
three analyses. 

Getting radiation forces right is notoriously difficult. 
Even for Cassini, whose construction is well known, the 
predicted and measured values differ by 50% [3]. How- 
ever, the total force can be no larger than the sum of 
the possible components, though it can easily be less. 
Therefore the main job is to show that enough force is 
available; any lesser result is easily explained. 



Figure 2 largely follows Figure 1 chronologically. Pio- 
neer 10 was at 40 AU in 1987.) The authors divide the 
history into three intervals. Interval I is January 1987 to 
July of 1990, interval II from July of 1990 to July of 1992, 
and interval III is from July of 1992 to the June of 1998. 
The authors make this distinction by looking at the spin 
rate of the craft (see Figure 3). In intervals I and III it 
was decreasing smoothly, but in interval II it decreased 
quickly and irregularly. They therefore consider the data 
from interval II to be less reliable than intervals I and III, 
since whatever affected the spin in interval II (probably 
gas leaks) may also have affected the acceleration. 



II. THE ANOMALOUS ACCELERATION 

As the Pioneer spacecraft recede from the sun, solar 
forces decrease and only gravitational forces, and an oc- 
casional maneuver, should affect the trajectory of the 
spacecraft. Anderson, et al. noticed that a small ad- 
ditional acceleration needed to be added to make the 
measured data and computations match. This is the 
anomalous acceleration, which started to become notice- 
able about 5 AU from the sun. It was roughly the same 
for Pioneer 10 and 11, as shown in Figure 1. 
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FIG. 1: Unmodelled acceleration as a function of distance 
from the sun, by Anderson et al. [3]. 



Additional constraints come from further study of Pi- 
oneer 10, since the data are higher quality and the data 
span is long enough to provide significant constraints due 
to the radioactive decay of the heat sources. Figure 2, re- 
produced from [4], shows the measured acceleration 1987 
to 1998. (Although they have different horizontal axes. 
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FIG. 2: Unmodelled acceleration and an empirical fit from 

Turyshev[4]. 



More recent analyses have refined these results some- 
what, though the main conclusions remain unchanged. 
Three different analyses have been reported in the litera- 
ture. SIGMA and CHASMP are two difli'crcnt trajectory 
modelling programs each with many possible analysis op- 
tions. We use the best Weighed Least Squares (WLS) 
results from each program, from [3]. Markwardt [2] wrote 
an new program with the explicit goal of an independent 
re- analysis. 

Table I shows the most recent results from [3], which 
fits a constant, independent acceleration in each inter- 
val. Table II shows the results of Markwardt's re-analysis 
which fits a constant plus a linear term to the data from 
1987-1994. His best solution is 

a{t) = -8.13 • lO-^cm/sec^ + 3.7 • 10~"t cm/sec^ 

where t is the time in seconds since the beginning of 
1987. Accelerations are in units of 10^* cm s^^. For 
convenience, we show the amount of directed power, in 
watts, that would be needed to account for each acceler- 
ation, assuming the 241 kg estimate of spacecraft mass 
from [3]. 

Note that each program claims very small formal er- 
rors, but the programs differ from each other by far 
greater amounts. Therefore the errors are probably sys- 



TABLE I: Summary of results from Anderson, et a/. [3] 



Interval 


SIGMA 
accel. 


equiv. 

watts 


CHASMP 
accel. 


equiv 

watts 


Jan 87- Jul 90 


8.00 ± 0.01 


57.8 


8.25 ± 0.03 


59.6 


Jul 92 - Jul 98 


7.84 ± 0.01 


56.7 


7.91 ± 0.01 


57.2 



TABLE II: Summary of results from Markwardt [2] 


date 


accel. 


equiv. 
watts 


Jan 87 - Mar 94 
(all data, constant ace.) 


7.70 ± 0.02 


55.7 


Jan 87 - Jul 90 
(constant ace.) 


7.98 ± 0.02 


57.7 


Jan 87 
(from linear fit) 


8.13 ± 0.02 


58.8 


Jan 87 - Jul 90 
(from linear fit) 


7.93 ± 0.02 


57.3 


Jul 92 - Jul 98 
(extrapolated from linear fit) 


7.14 


51.6 



tematic, not random, and the differences between the 
programs are better estimates of the real uncertainties. 
If we create a meta-analysis by averaging over the 3 anal- 
yses (rather dubious, but the best we can do) we get 
an acceleration in interval I of 8.08 ± 0.12 (58.3 ± 0.87 
watts). Agreement here is good, so this number can be 
regarded as fairly secure. The variation with time is less 
clear, with SIGMA and CHASMP showing a 2.00% and 
4.12% decrease between interval I and interval III, and 
Markwardt finding a linear trend that predicts a 10.6% 
decrease in this interval. This gives a meta-prediction of 
(5.55 ± 3.63)% decrease from interval I to III. 

Although much less publicized, there are other unmod- 
elled forces acting on the craft as well. In the absence of 
external forces and/or spacecraft structure changes, the 
spin rate should not change. It does change, though, as 
shown in Figures 3 and 4, from[3]. Note that Pioneer 
10 is spinning down at a rate proportional to the bus 
power (in intervals I and III), and Pioneer 11 is spin- 
ning up, except at manuevers. From the viewpoint of 
fundamental physics, unexplained torques are as inter- 
esting as unexplained forces, and the two are of compa- 
rable size (13 watt-meters and 57 watts for Pioneer 10). 
Nonetheless, [3] assumes that the spin changes are caused 
by spacecraft systematics, but tries to show the acceler- 
ation changes are not. This distinction is driven partly 
by the data (the spin rates change at boundaries defined 
by spacecraft events, the acceleration does not) but also 
by a lack of any remotely plausible alternative. There 
are many (and interesting) theories that could cause ac- 
celeration (modified gravity, dark matter, and so on) but 
there are few proposed theories that could cause anoma- 
lous spins. 
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FIG. 3: Spin changes in Pioneer 10 [1]. The vertical lines 
indicate the times of maneuvers. 
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FIG. 4: Spin changes in Pioneer 11 [1]. The vertical lines 
indicate the times of maneuvers. 



III. PREVIOUS WORK 

For the convenience of the reader, section III A consists 
of direct quotes from [3] , covering the relevant details of 
the Pioneer spacecraft, and Figure 5, from [5]. Many 
other paper[5] and web[6, 7] descriptions are available. 
In section III C we summarize the existing literature on 
the hypothesis that non-isotropic radiation is responsible 
for the unmodellcd acceleration. 



A. General description of the Pioneer spacecraft, 
from [3] 

The main equipment compartment is 36 cm deep. 
The hexagonal flat top and bottom have 71 cm long 
sides. Most of the scientific instruments' electronic units 




SEPARATtON 

INFiiAi!Ei> ■ 

' fCil/MUNCm . 

CHARGED PART! CLE . 

iNStRUMEt^T 

(CH !CAO 0/SIMPSON} 

MAGNETOMfTER . 
fDEPLOVeD) fJPt/SMfTH) 

TKAf PED RAOiAtlON 
/ .' DETECTQR <UC5G/Fll!.!LISS 

arSAVIt^ET 
PHOTOMfTJft 

LOW ^AIN 
/.(^T^NA 



STA 

-17. ?0 

ASTfRO!&/KiiEtEOROiU 
DSTFCTOE (GE^OffERMAN) 
>CGSPACEOWFT 
GEKSERTUBg 
TElE$COf>e 
(lOWAA'AN AUENS 

IMAGING 

PHOTO?Ot«!lfcffTH! 

tU/ARIZONA/G£Hi!ELS) 



FIG. 5; Reproduction of Figure 3.1-2 from [5]. A few lines 
were removed for clarity, and the main equipment compart- 
ment is shaded in. 



and internally-mounted sensors are in an instrument bay 
( "squashed" hexagon) mounted on one side of the central 
hexagon. 

At present only about 65 W of power is available to Pi- 
oneer 10 [8]. Therefore, all the instruments are no longer 
able to operate simultaneously. But the power subsystem 
continues to provide sufficient power to support the cur- 
rent spacecraft load: transmitter, receiver, command and 
data handling, and the Geiger Tube Telescope (GTT) sci- 
ence instrument. 

The sunward side of the spacecraft is the back, and 
the anti-sunward side, in the direction of motion, is the 
front [9]. 



B. Gas leaks 

Gas leaks are always a prime suspect when unmodelled 
spacecraft accelerations are found. As the authors them- 
selves say "Although this effect is largely unpredictable, 
many spacecraft have experienced gas leaks producing 
accelerations on the order of 10~^ cm/s^"[3]. Further- 



more the authors think that gas leaks are a significant 
part of the spin behavior. Why then do they think that 
gas leaks are not the source of the acceleration? They 
present four arguments: 

• The effect seems constant over long periods of time 
(many years). 

• The acceleration does not change as a result of 
thruster activity, as many gas leaks do (as valves 
seat/unseat.) 

• The effect is roughly the same on two spacecraft. 
Pioneer 10 and 11. 

• A force big enough to cause the acceleration would 
cause bigger spin changes than are observed unless 
it was directed along the spin axis. 

In rebuttal, there are many possible sources of gas leaks, 
not all of which are variable or affected by thruster ac- 
tivity. (The same authors [3] speculate that a gas leak 
causes the spin- up of Pioneer 1 1 , which is also constant 
and unaffected by maneuvers through 4 years.) Further- 
more, the two spacecraft were intended to be identical, so 
an identical artifact such as a gas leak would not be sur- 
prising, and it might be aligned with the axis. In short, 
it would require an unusual gas leak, duplicated on each 
spacecraft, to cause the observed effect, but it is certainly 
allowed by physics. 

In [3], the error budget for gas leaks is set as follows: 
First, take the biggest uncommanded spin-rate change, 
assume it was caused by gas leaks, assume the leak was 
at the spin thrusters, and then increase it a little. Thus 
they are setting the budget to the biggest known leak 
on this particular spacecraft. This is hardly a rigorous 
method for estimating the maximum possible size of an 
unknown leak, since there could be more than one leak, 
and locations other than the thrusters require a bigger 
leak for the same spin change. Furthermore, as the au- 
thors note, other spacecraft are known to have had larger 
leaks. Clearly at least some of the authors of [3] are not 
convinced by their own argument since they still suspect 
gas leaks as the cause of the unmodelled acceleration ([3], 
Section XII). 



C. Non-isotropic radiation - previous work 

Murphy suggests that the anomalous acceleration seen 
in the Pioneer 10/11 spacecraft can be, "explained, 
at least in part, by non-isotropic radiative cooling of 
the spacecraft." [10], The main idea is that heat from 
the main and instrument compartments would radiate 
through the cooling louvers on the front of the craft. An- 
derson, et al. argue in reply [11] that over the data span 
in question the louver doors were already closed (if the 
doors were open then the effect would surely be signif- 
icant). They conclude "the contribution of the thermal 
radiation to the Pioneer anomalous acceleration should 
be small." They also argue that the spacecraft power 
is decreasing, but the unmodelled acceleration is not. 



SchefFcr[12] points out that the front of the spacecraft 
has a much higher emissivity than typical thermal blan- 
kets (even with the louvers closed), and therefore the 
majority of the heat will radiate from the front in any 
case. Anderson, et al. [13] dispute this, based on the 
emissivity data in [5] , which assigns a high emissivity to 
the thermal blanket. (This data is true but misleading - 
it specifies the emissivity of the outer layer of the blanket. 
This is very different from the emissivity of the blanket 
as a whole, called the effective emissivity, which is quite 
low. The next section has a more in-depth discussion of 
this point.) 

Katz[14] proposes that at least part of the acceleration 
is generated by radiation from the RTGs reflecting off the 
back of the antenna. Anderson et al. in [15] argue that 
this effect must be small since the antenna is end-on to 
the RTGs, and hence gets very little illumination. 

Slusher (as credited by Anderson) proposed that the 
forward and backward surfaces of the RTGs may emit 
non-equally. Anderson et al. conclude there is no credible 
mechanism to explain the large difference in surfaces that 
would be required if this was to explain the whole effect. 



IV. DISCUSSION 

We consider asymmetrical radiation from 5 sources - 
the RTGs themselves, the two spacecraft compartments, 
RTG radiation reflected from the antenna, the Radioiso- 
tope Heater Units (RHUs) on the spacecraft, and radi- 
ation from the feed that misses the antenna. We also 
consider one modelling error, a mis-estimation of the re- 
flectivity of the antenna to solar radiation. 

Consider thermal radiation from the spacecraft body 
with the louvers closed, as they have been since 9 AU. 
An extremely simple argument shows that the electri- 
cal power dissipated in the main spacecraft compartment 
must result in a significant amount of thrust. The Pi- 
oneer antenna points roughly at the sun, and the in- 
strument compartment is directly behind the antenna. 
Since the antenna blocks radiation in the sunward di- 
rection, the waste heat must be preferentially rejected 
anti-sunward. Referring to Figure 5, a good scale model 
is a 60 watt bulb about 4 cm behind a 25 cm diameter 
pie dish. The dish casts a huge shadow in the sunward 
direction, resulting in an average anti-sunward thrust. 

However, the efficiency of conversion of heat to thrust 
is higher than this simple argument indicates. From [5], 
"The Pioneer F/G thermal control concept consists of 
an insulated equipment compartment with passively con- 
trolled heat rejection via an aft [25] mounted louver sys- 
tem." Since even a closed louver is a much better radia- 
tor than thermal insulation, most of the radiation occurs 
from the front. It's as simple as that! 

Instrument heat may also contribute to thrust, but 
possibly with less efficiency. This is because the instru- 
ments could possibly radiate at right angles to the spin 
axis through their observation ports, which are not cov- 



ered with thermal blankets. Furthermore, the science 
compartment is much closer to the edge of the dish than 
the main compartment, so the dish will shadow much 
less of any thermal radiation generated by the science 
instruments. 

We estimate the efficiency using the spacecraft con- 
struction. Assuming a uniform internal temperature, the 
power emitted from each surface is proportional to the 
area times the effective emissivity of the surface. The 
front and back of the central equipment compartment 
have about 1.3 m^ area, and the sides about 1.5 w? to- 
tal. The sides and the rear of the compartment are cov- 
ered with multi-layer insulation(MLI)[5]. When calcu- 
lating radiation from multi-layer insulation, the correct 
value to use is the "effective" emissivity, Ce//, which ac- 
counts for the lower temperature of the outer layer[16]. 
(Anderson[13] points out that the outer layer of the MLI 
has an emissivity of 0.70 according to [5]. This is not 
a contradiction because the outer layer of the MLI is 
much colder than the interior - that's how MLI works.) 
From [16], the multilayer insulation from on Pioneer 10 
has an effective emissivity of 0.007 to 0.01 (see Figure 
6). Assuming a value of 0.0085, and a 1998 internal tem- 
perature of 241 K [17], the main compartment will lose 
about 4 watts total through the MLI on the sides and 
back (Even this may be an over-estimate. Two of the 
sides are facing 1 Kw IR sources just 2 meters away, and 
may even conduct heat into the compartment.). Allow- 
ing a few watts for conduction losses through wires and 
struts, perhaps 10% of the power (about 6 watts) goes 
through the back, 10% through the sides, and the remain- 
ing 80% through the front. The back radiation will have a 
near zero efficiency (it squirts out from between the dish 
and the compartment at right angles to the flight path). 
Radiation from the side should be about 10% efficient, 
assuming Lambcrtian radiation and a 45 degree obstruc- 
tion by the dish. Radiation from the front will be about 
66% efficient, again assuming Lambertian emission. The 
overall efficiency of main bus radiation could therefore be 
as high as 54%. 

Is it reasonable for the front of the main compartment 
to radiate the 47 watts or so this requires? At an aver- 
age temperature of 241 K, and assuming a flat surface, 
this would require an average emissivity of 0.19. From a 
picture of the Pioneer 10 replica in the National Air and 
Space Museum [18], the front of the spacecraft is rather 
complex, with considerable surface area (such as a rather 
large cylinder that connects to the booster) and a variety 
of surface flnishes. There are also some fairly large instru- 
ments on the front of the spacecraft, such as the plasma 
analyzer [5] . Although the louver blades themselves have 
a low emissivity of 0.04 [5], a composite emissivity of 0.19 
seems reasonable. 

The main conclusion seems quite robust. Multi-layer 
insulation is specifically designed to reduce heat losses, 
whereas the louvers have at most one layer of obstruction 
even when closed, and by definition are riddled with dis- 
continuities, which are a major source of heat leaks [16]. 
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For the Pioneer spacecraft, P — 8 watts and A = 13 cm, 
so the torque is about 0.165 watt-meters. This is about a 
1/100 of the total observed torque and can be neglected. 



B. Radiation from the RHUs 

From the diagrams in [5], 11 1-watt (in 1972) radioiso- 
tope heater units are mounted to external components 
(thrusters and the sun sensor) to keep them sufficiently 
warm. The diagram is not very specific, but it appears 
that 10 of the units, those mounted to thrusters and the 
sun sensor, are behind the edge of the main dish. If we 
assume these radiate isotropically into the hemisphere 
behind the antenna, then they contribute the equivalent 
of 4 watts of directed force in 1998. The remaining RHU 
is at the magnetometer and would not appear to con- 
tribute net thrust. RHU thermal radiation will decrease 
with a half life of 88 years. 



C. Asymmetrical radiation from the RTGs 



FIG. 6: Figure from [16]. Pioneer F became Pioneer 10. 



The lowest emissivity material on the front, 0.04, has at 
least 4 times the highest quoted emissivity of the sides 
and back. Surely, therefore, a majority of the heat will 
be radiated from the front of the spacecraft. 



A. Feed pattern of the radio beam 

An ideal radio feed antenna would illuminate its dish 
uniformly, with no wasted energy missing the dish. How- 
ever, the feed is physically small and cannot create such a 
sharp edged distribution, so some radiation always spills 
over the edge. Since dish area is wasted if not fully il- 
luminated, an optimum feed (for transmission) normally 
allows about 10% of the total power to miss the dish. 
This power is converted to sunward thrust with an effi- 
ciency of 0.7 since it is directed roughly 45 degree angle 
to the spin axis. The rest of the energy hits the antenna 
and is reflected sunward. If epEED is the fraction of the 
energy that misses the antenna, and the transmitter is 8 
watts, then the net thrust towards the sun is 

(8 w)(eFEED • 0.7 - (1 - epEEo)) 
This is negative as expected, since most of the radiation 
is sunward. 

As a side note, the radio beam is circularly polarized 
and thus carries angular momentum away from the space- 
craft. A circularly polarized beam of power P and wave- 
length A will impart a torque T of 

T^p^n 

he 



The RTGs might contribute to the acceleration by ra- 
diating more to the front of the spacecraft than the rear. 
This might be caused by differing solar wind and/or dust 
enviroments, as proposed by Slusher. In [3] this is ana- 
lyzed and estimated to be a small effect. 

However, a more likely cause has recently been pro- 
posed by Herbert [19], who notes that solar UV radiation 
can bleach coatings such as that used on the RTGs, and 
hence make the sunward side a slightly worse radiator 
than the unexposed anti-sunward side. There is at least 
some experimental evidence to indicate this is plausible. 
Two commonly used coatings, Z93 and YB71, broadly 
similar to those on the Pioneer RTGs (metal oxide pig- 
ment and a silicate binder), were tested on the NASA 
Long Duration Exposure Facility (LDEF) which was left 
in orbit for several years and then retrieved by the space 
shuttle. The IR reflectance/absorbance characteristics of 
the coatings were measured after the exposure to orbital 
conditions. The Z93 coating was essentially unchanged 
as an IR emitter (see Figure 49 of [20]), but the YB71 
coating was a worse emitter, by up to few percent depend- 
ing on the wavelength (Figure 57 of [20]). Weighting this 
curve by the blackbody spectrum of a 440K radiator gives 
about a 3.2% degradation in IR emission. This evidence 
is certainly not proof that this effect has occurcd, since 
the exact coatings and environments differ (for example, 
the earth orbit encounters more charged particles, atomic 
oxygen, and contamination, in addition to the solar UV), 
but indicates that it is a plausible suspect. 

In [3], RTG asymmetry is found to give about 6 watts 
of thrust for a 1% asymmetry, which they take as the 
maximum plausible value. From the data above, this 
maximum is too low, since 3.2% has been observed on 
similar surfaces. In this paper we will assume that 1% 
bleaching occurred, resulting in 6 watts of thrust, but up 



to 3% is certainly plausible in the absence of additional 
data. 

Asymmetrical RTG radiation (to one side, not fore and 
aft) could also be the cause of the slow but constant spin- 
up observed on Pioneer 11. This would explain why the 
spin-up is almost constant in value and unaffected by 
manuevers. 



D. Revisiting RTG reflection 

Some of the waste heat from the RTGs will reflect from 
the back of the high gain antenna and be converted to 
thrust, as proposed by Katz[14]. Anderson et al.[15] ar- 
gue that at most 30 watts of radiation hit the antenna, 
and hence RTG reflection cannot account for the whole 
acceleration, which is true. Similarly, Slabinski[21], in an 
unpublished analysis from 1998, concluded that roughly 
28 watts of radiation hits the antenna, and hence the 
whole effect could not be explained. However, it is clear 
the effect is real, and can provide a significant fraction 
of the observed anomaly. Only the exact amount is in 
question. 

The RTGs are not on-axis as viewed from the antenna. 
From Figure 5, we see that the centerline of the RTGs is 
behind the center of the antenna. Measurements from 
this diagram indicate this distance is about 23.8 cm. 
(Slabinski[21] independently estimated 26 cm for this dis- 
tance.) Another figure (not included here) from [5] shows 
the far end of the RTGs is 120.5 inches (or 3.06 meters) 
from the centerline. The near end of the RTGs will then 
be about 60 cm further in, or at about 2.46 meters from 
the center. The antenna extends 1.37 meters from the 
center, so the rim of the antenna is 69.8 cm off axis and 
1.09 meters away radially. Thus the edge of the antenna, 
where the illumination is by far the brightest, views the 
inner RTG at an 32.6 degree angle. This is far from on- 
axis. 

The fins of the RTGs radiate symmetrically, and all are 
visible from the antenna, so the center of this illumina- 
tion will be 23.8 cm behind the antenna. The cylindrical 
center of the RTG is about 8.4 cm in radius [22] so this 
illumination will come from at about 15.4 cm behind the 
antenna. The fins have more area than the cylinder, so 
for this calculation we take a rough weighted average and 
assume a cylindrical Lambertian source 20 cm behind the 
antenna. We assume the inner RTG is centered 2.66 me- 
ters from the center, and the outer RTG 2.91 meters. 

The area blocked by the antennas is shown in Figure 

7 in spherical coordinates. Numerical integration of the 
two areas shows about 12 watts for the near RTG and 

8 watts for the far one if the total RTG power is 2000 
watts. This does not include radiation from the endcaps 
or supporting rods. 

Combining the analyses, we conclude that at least 20 
watts, but no more than 30 watts, of radiation hits the 
antenna. In this paper, we will use 25 watts as the basis 
for further analysis. This energy is turned into thrust by 




FIG. 7; Antenna size in spherical coordinates from RTGs. 
Radial axis is angle from the centerline in radians; other axis 
is angle around this line with the magnetometer defined as 



two effects. First, the antenna shadows radiation which 
would otherwise go forward. An angle in the middle 
of the antenna is about 17 degrees forward; this corre- 
sponds to an efficiency of 0.3 (the true efficiency is prob- 
ably higher since the edge is both at a greater angle and 
more brightly illuminated.) Next, the energy that hits 
the antenna must go somewhere. Some will be absorbed 
and re-radiated; some will bounce into space, and some 
will bounce and hit the instrument compartment, and be 
reflected or re-radiated from there. A detailed account- 
ing seems difficult, but an overall efficiency of 0.6 to 0.9 
seems reasonable (0.3 for shadowing and 0.3 to 0.6 for 
reflection and re-emission). 



E. Total of all eff'ects 



Here we sum the maximum value of all the effects as 
of 1998. The total is more than enough to account for 
the acceleration, giving us the freedom to reduce some of 
the efficiencies if needed to fit the data. 



TABLE III: Available thrust from different sources as of 1998 



Source of effect 


Total Power 


Effic. 


Thrust 


Decay 


Rad from RHUs 


8 


0.5 


4 


0.78%/year 


Antenna shadow 


25 


0.3 


7.5 


0.68%/year 


Antenna radiate 


25 


0.6 


15 


0.68%/year 


RTG asymm. 


2000 


0.009 


18 


0.68%/year 


Feed pattern 


0.8 


0.7 


0.6 


0% 


Radio beam 


7.2 


-1 


-7.2 


0% 


Rad., main bus 


59 


0.54 


32 


see text 


Rad., instr. 


1 


0.1 


0.1 


see text 


Total 






70.0 





F. Antenna solar reflectivity 

In this section wc argue that a mismodelled solar reflee- 
tion might account for the sudden onset of the anomalous 
force shown in figure 1. This argument is offered only as 
a possible explanation of the onset and initial decrease of 
the anomalous acceleration; it is not relevant to the ex- 
istence, magnitude, or source of the acceleration at later 
times since past about 30 AU the contribution from the 
solar radiation is negligable. 

First, we show there is surely a possibility of error in 
these coefficients since the numbers for the two space- 
craft disagree. We start with the data from Anderson, 
averaging the SIGMA and CHASMP values. We assume, 
following Anderson (section VII-B), that the trajectory 
was fit correctly but the mass used in the calculation 
was incorrect. We correct the fitted values using the best 
available estimates for spacecraft mass, keeping the ac- 
celeration (and hence trajectory) the same. We would 



TABLE IV: Solar reflectivity from Anderson, et al. 



Spacecraft 



Pioneer 10 



Pioneer 11 



fitted K. 



1.73 



1.83 



mass used 
in fit 



251.8 



239.7 



actual 
mass 



241 



232 



resulting 
true /C 



1.66 



1.77 



expect an nearly identical value of /C for both spacecraft 

- they were the same size and painted with the same paint 

- but we observe /C = 1.66 for one and K, ~ 1.77 for the 
other, a 6.6% difference. One possible explanation is that 
the two spacecraft had different amounts of thermal radi- 
ation thrust, and the fitting procedure used an adjusted 
value of /C to fit the observed trajectory. 

Is it possible that the trajectory is right, but /C is 
wrong? Anderson et al.[13\ claim that any appreciable 
error in this value would have resulted in navigation er- 
rors, but the 6.6% difference between the Pioneers (surely 
not physically present) was easily absorbed into the fit 
(perhaps into the velocity increments at maneuvers, for 
example). This would certainly not be the first time that 
an excellent fit to the data was obtained with the wrong 
explanation. 

The analysis in Anderson [3] reinforces this point. They 
tried to separate a constant anomalous force from the 
1/r^ solar reflectivity, but found they were tightly corre- 
lated. On Ulysses, for example, the correlation between 
these two parameters was 0.888, so 90% of the change in 
one parameter could be explained by a spurious change 
in the other. If the hypothesis of this article is correct, 
the two parameters will be even harder to separate, since 
for Pioneer the both the radiation contribution and the 
solar reflection are decreasing as the spacecraft recedes 
from the sun (the total power is decreasing, and the effi- 
ciency will decrease as well as the louvers close). 

In the scenario of this paper, the acceleration has ex- 
isted all along, and might even have been stronger closer 



to the sun. When Pioneer was closer to the sun, though, 
the fitting programs absorbed the extra acceleration by 
adjusting the value of /C and perhaps other parameters 
such as the delta-v of maneuvers. As Pioneer receded 
from the sun, and maneuvers became less frequent, ad- 
justments to these parameters could no longer fit the tra- 
jectory properly. At this point the the anomalous accel- 
eration "appears" . This argument is not specific to radi- 
ation induced acceleration - any small radial acceleration 
can be compensated for by adjusting the value of /C. 

In this paper, we model the effect of any error in /C by 
introducing a fictitious force, whose value is simply the 
solar force on the spacecraft times the error in /C. 



V. COMPARISON WITH EXPERIMENT 

How well does this explanation account for the accel- 
eration? The explanation has 6 adjustable parameters. 

• ^RHU I the fraction of RHU heat converted to thrust 

• ^RTGi the fraction of RTG heat converted to thrust. 
Includes both direct asymmetry and reflection from 
the antenna. 

• ^FEED, the fraction of RF power that misses the 
antenna 

• ^INST, the fraction of instrument heat that is con- 
verted to thrust 

• ^Bus I the fraction of main compartment heat that 
is converted to thrust 

• K SOLAR,, the amount by which the solar reflection 
constant is underestimated. This cannot exceed 
about 0.2 since the true value can be no more than 
2.0. 

In theory all parameters are separable since they decay 
at different rates. In practice there are many similar 
solutions that cannot be distinguished by the existing 
data. 

We compute the net thrust as follows: let d be the date 
in years. The total electrical power, in watts, is modelled 
as 

E{d) == 68 + 2.6- (1998.5 -d) 
The RHU power, in watts, is 

RHU{d) = 10.0 • 2-('^-i972)/88 

The RTG heat dissipation, in watts, is 

RTG{d) ^ 2580 ■ 2-(''-i9^2)/88 _ ^^^^ 

We assume the distance from the sun, measured in AU, 
increases hnearly from 20 AU in 1980 to 78.5 AU in 2001: 

r{d) =2Q+{d- 1980)/21 • (78.5 - 20) 

The power incident upon the antenna, in watts, is 

SOLAR{d) = 7r(1.37 mffQ/r^{d) 

where /q = 1367 W/m (AU)^ is the "solar radiation 
constant" at 1 AU. We use the expression from section 
IV A for the radio thrust. 

The power dissipated in the instrument compartment, 
INST{d), is given in Table V from [23]. 



TABLE V; Instrument power 1987-2001. IPP = Imaging Pho- 
topolarimeter, TRD = Trapped Radiation Detector, PA = 
Plasma Analyzer 



Dates 


Watts 


Notes 


Jan 87 - Oct 93 


11.6 


IPP off Oct 93 


Oct 93 - Nov 93 


8.1 


TRD off Nov 93 


Nov 93 - Sep 95 


5.3 


PA off Sep 95 


Sep 95 - present 


0.8 


Only Geiger active 



The other units that were turned off during this period 
(the Program Storage and Execution unit, and the Du- 
ration and Steering Logic) did not affect the instrument 
heat since they simply substituted one heat source in the 
main compartment for another. 

The electrical power that does not go into the instru- 
ments or the radio beam goes into the main compart- 
ment: 

BUS{d) = E{d) - INST{d) - 8.0 

We sum the individual sources, then convert to accel- 
eration by dividing by c, the speed of light, and m, the 
spacecraft mass (here 241 kg): 

accid) = \eRHU ■ RHU(d) + 

c ■ m 

CRTG ■ RTG{d) + 

(8 w)(eFEED • 0.7 - (1 - efeed)) + 

eiNST ■ INST(d) + 

CBUS ■ BUS{d) - 

KsoLAR * SOLAR{d)] 
To examine the fit, we use the plots from [3, 4], and try 
to fit them with our model. We make three fits. The first 
is a conservative fit, using only known and documented 
spacecraft characteristics. The second is the nominal fit, 
adding in effects such as RTG asymmetry that are plau- 
sible but not proven. The third is constructed to get 
the best possible fit to the data, but might be physically 
unrealistic. 

The conservative fit uses only known and documented 
spacecraft characteristics. These are that the front of the 
spacecraft is a better radiator than the sides [5, 16], and 
that the antenna will block and reflect some of the RTG 
radiation[3, 5, 14, 21]. A good fit is obtained with 

• ^RTG = 0.01. 25 watts hit the antenna, 30% block- 
age efficiency, and 50% reflection efficiency. 

• ^INST = 0.51. Instruments same as main bus for 
simplicity. 

• ^BUS = 0.51. About 80% the main bus heat goes 
out the front, with Lambertian efficiency. 

This model correctly predicts 58.6 watts in interval I, 
but predicts a decrease in interval III to 48.6 watts. This 
is a 17% decrease as opposed to the 3% measured in 
Anderson and 10.6% of Markwardt. This model does not 
explain the onset at 5 AU, and overpredicts the rate of 
decrease, but it shows that at most 20% of the effect can 



be due to new physics. At the very least, 80% of the effect 
can be accounted for by entirely conventional physics, 
based on known, documented, and measured spacecraft 
construction. 

The nominal fit adds radiation from the RHUs, asym- 
metrical radiation from the RTGs, feed spillover, and so- 
lar reflectance mis-modelling. These sources are all plau- 
sible but neither proven or disproven by any records or 
measurements found so far. The fit assigns the same ef- 
ficiency to main compartment heat and instrument heat. 
This avoids much of the need to look at spacecraft con- 
struction details and instrument history, since the accel- 
eration only depends on the total electrical power. 

The additional sources allow a better fit since RHU and 
RTG heat decays more slowly than electrical heat, feed 
spillover does not decay at all, and we can now model 
the onset of the anomalous acceleration at 5 AU. Once 
again, many parameter choices give similar results. We 
get a reasonable fit over the entire data span with the 
following coefficients: 

• ^RHU ~ 0.5, the RHUs radiate like point sources 
behind the antenna. 

• eRTG = 0.016. 0.3% RTG asymmetry, 30% block- 
age efficiency, and 50% reflection efficiency. 

• ^FEED — 0.1. 10% of the feed power misses the 
antenna 

• ^INST = 0.39. Instrument heat radiates as main 
bus heat for simplicity. 

• ^BUS = 0.39. About 60% of the main bus heat goes 
out the front, with Lambertian efficiency 

• KsoLAR — 0.2. Antenna reflection estimates are 
too low by 0.2. 

The fit to the data is shown in Figures 8 and 9. The 
agreement seems reasonable in both regimes. In partic- 
ular, the early anomalous acceleration between 15 and 
40 AU is fit well by this model. In Figure 8 two other 
models are shown, all assuming that a 1/r^ error of some 
sort (here solar constant mismodelling) is responsible for 
the onset. The middle trace assumes the acceleration is 
a pure exponential with an 88 year half life. This is the 
form for a model that assumes RTG radiation (direct or 
reflected) is asymmetric but spacecraft electrical heat is 
radiated isotropically. Between 15 and 40 AU this model 
underpredicts the observed decrease, where the nominal 
model flts much better. This strongly favors a model 
where radiation from the spacecraft bus is a major con- 
tributor to the anomalous acceleration. The lower trace 
is a constant acceleration plus an error that scales as 
l/r^. This shows that if the acceleration is indeed con- 
stant at large distances, a different explanation for the 
onset is required. 

The flt from 1987 to 1998, shown in Figure 9, also 
looks reasonable. We compare this model to the consen- 
sus of the most recent analyses [3] [2]. Using the param- 
eters above, the average sunward thrust is 58.0 watts in 
interval I, and 50.2 watts in interval III. We can adjust 
the the parameters to get the correct overall average, or 
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FIG. 8: Data from Figure 1 (error bars), model prediction 
from this paper (solid line), pure 88 year half life plus solar 
constant error (middle line), and constant acceleration plus 
solar constant error(lowest line) 
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FIG. 9; Figure from [4], with fitted data added. The dotted 
line is Turyshev's empirical fit; the solid line is the model 
hypothesized in this paper. 



the right acceleration in interval I, but in either case we 
would expect to see a 13.2% decrease from interval I to 
III, where only a 5.6% decrease is observed. The 7.6% 
discrepancy is about 2 standard deviations out. Taken 
at face value, this makes it unlikely at about the 2% 
level that this hypothesis alone accounts for all the mea- 
sured result. However, the re-analysis by JVIarkwardt [2] 
has concluded that the data does not rule out a slowly 
decreasing force, at least if the the decrease is an ex- 
ponential with a half life of more than 50 years. (This 
corresponds to a 9% decrease in the 6.75 year span be- 
tween the midpoints of intervals I and III) . The decrease 
here is not strictly exponential, but is close in in size and 



shape to that explicitly allowed by Markwardt. 

JVIore speculatively, an even better fit to the accelera- 
tion data can be obtained by assigning different efficien- 
cies to instrument heat and main compartment heat. For 
example, 

• ^RHU ~ 0.5, the RHUs radiate like point sources 
behind the antenna. 

• ^RTG = 0.01425. 0.3% RTG asymmetry, 30% 
blockage efficiency, and 60% reflection efficiency. 

• ^FEED = 0.1. 10% of the feed power misses the 
antenna 

• ^INST — 0.40. About half the main bus heat radi- 
ated forwards with Lambertian efficiency. 

• ^Bus = 0.10. Instruments radiate mostly to the 
side 

• KsoLAR — 0.2. Antenna reflection estimates are 
too low by 0.2 . 

gives a better fit, with only a 4.9% discrepancy (10.5% 
predicted versus 5.6% measured) on the I-III decline and 
a roughly equivalent fit at earlier times. This is only 
about 1.3 standard deviations from the consensus model, 
and an almost perfect fit for Markwardt. However, fig- 
uring the maximal reasonable difference between instru- 
ment efficiency and and main compartment efficiency is 
difficult [13]. One the one hand the two compartments 
are separate, the instrument bay is closer to the edge of 
the antenna, and it has has side facing ports that extend 
though the thermal blankets. On the other hand the two 
compartments are radiatively and conductively coupled. 
Without a much more detailed analysis it's very hard to 
determine the maximum plausible difference in efficien- 
cies. 

Finally, asymmetric radiation offers a parsimonious ex- 
planation for both the anomalous acceleration and the 
anomalous torque. Anderson et al. note that the Pioneer 
10 spin-down torque is almost perfectly correlated with 
the main bus power. Radiation from the front of the 
craft, as proposed here, explains this. The needed emis- 
sion geometry is numerically plausible - in 1986, there 
were 97 watts available, and 13 w-meters of torque mea- 
sured. Assuming the radiation is emitted 50 cm from the 
axis (the louver location) , if the radiation was canted at 
an average angle of 15.5 degrees from the normal to the 
surface, it could provide the observed torque. Such an 
angle would decrease the conversion of power into thrust 
by only 4%, leaving that argument intact. The louvers, 
covering the front surface and all canted to one side when 
closed, provide a natural explanation for the asymme- 
try required. One obvious objection to this explanation 
is that it predicts Pioneer 11 should be spinning down 
as well, instead of the spin-up that is actually observed. 
This is not a serious problem since the unknown spin-up 
mechanism, possibly gas leaks or RTG asymmetry, can 
easily overpower the small torque induced by main bus 
radiation. 

In any case, the proposed explanation, by accounting 
for the bulk of the effect, makes it more likely that con- 
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ventional physics can account for the entire unmodelled 
acceleration. Conventional explanations for the remain- 
ing discrepancy include other unmodelled effects such as 
gas leaks, inaccuracies in the simple thermal model, or 
the effects of a complex fitting procedure applied to noisy 
data. 



VI. CONCLUSIONS AND FUTURE WORKS 

No new physics is needed to explain the behavior of the 
Pioneer spacecraft. Either gas leaks or thermal radiation, 
or a combination of the two, could explain both the linear 
and angular accelerations that are measured. 

A strong thermal effect is certainly present, based only 
on the construction of the Pioneers. Estimates show it 
can account for the magnitude of the unmodelled accel- 
eration to within the errors, but ovcrpredicts the rate of 
change. The antenna shadowing of main compartment 
radiation and the radiation from the RTGs falling on the 
antenna seem particularly robust sources of acceleration 
since they are only based on geometry. These effects 
alone account for more than half the acceleration. The 
other sources - RHU radiation, differential RTG radia- 
tion, and differential emissivity - depend more on con- 
struction details, but all seem plausible. 

This explanation also explains some other puzzles: the 
values of acceleration of Pioneer 10 and 11 would be 
expected to be similar, but not identical, as observed. 
The acceleration and the observed torque (on Pioneer 10) 
share a common origin, and the torque is proportional to 
the main bus heat, as observed. Other spacecraft, built 
along the same general principles, would be expected to 
show a similar effect, but planets and other large bodies 
would not, as is observed. 

The hypothesis here predicts an eventual, unambigu- 
ous decrease in the anomalous acceleration. If the accel- 
eration remains constant, on the other hand, the hypoth- 
esis will be refuted. Extending the analysis of Markwardt 
to the whole Pioneer data span would be useful, since it 
it currently stops at 1994 and it directly includes the 
possibility of a non-constant acceleration. Extending the 
analysis of Anderson by including post 1998 data would 
be helpful as well. 



If Pioneer 10 remains operational, additional data may 
allow us to improve our understanding of the unmodelled 
acceleration. The difference between constant accelera- 
tion and the decrease predicted by the hypothesis of this 
paper grows quadratically with time. Since the beginning 
of data in 1987, by 2002 the two solutions differ by 4.4 
cm/sec, or a doppler shift of 0.58 Hz. Thus a single good 
2002 measurement could tell the two hypotheses apart. 
Unfortunately the signal is now very weak, to the point 
where the standard JPL receivers have trouble locking 
onto the signal [7]. Careful recording of the return signal 
might probably work, though, with the frequency recov- 
ered through long averaging. Bigger telescopes such as 
Arecibo, the VLA, or Greenbank, might conceivably be 
pressed into service as well. 

More detailed modeling, using the Pioneer materials, 
construction details, and history, could provide a much 
better estimate of the magnitude of this effect. A suit- 
ably detailed thermal model, measured in a cold vac- 
uum chamber, would provide the strongest evidence for 
or against this hypothesis. 

Longer term, other proposed experiments such as 
LISA [24] are designed specifically to reduce the system- 
atics that bedevil retrospective analyses like Pioneer. 
(LISA is expected to be about 10'"' times better in this 
respect.) If the anomalous acceleration is not detected in 
these more precise experiments, then almost surely the 
unmodelled acceleration of Pioneer 10 is caused by over- 
looked prosaic sources such as those proposed here. 
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